Fractional periodicity and magnetism of extended quantum rings 
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The magnetic properties and nature of the persistent current in small flux-penetrated t — t' — U 
rings are investigated. An effective rigid-rotator description is formulated for this system, which 
coincides with a transition to a ferromagnetic state in the model. The criteria for the onset of effective 
rigid rotation is given. The model is used to understand continuum model ground-state solutions 
for a 2D few-particle hard-wall quantum dot, where ferromagnetic solutions are found even without 
the Zeeman coupling to spin. After the onset of effective rigid rotation, a 97-98% correspondence 
can be determined between the lattice model and continuum model eigenstate results. 

PACS numbers: 71.10.Fd, 73.23.Ra, 73.21.La, 74.25.Ha 



I. INTRODUCTION 

Electron traps in semi-conductor materials^ are 
nanosystems that are of intense current interest. Ex- 
amples include quantum rings and dots, which, due to 
their atom- like features, have shown immense potential 
for technological application. Recently, it has been a 
goal to acquire a more fundamental understanding of 
these systems. For quantum rings, this interest has 
been partly motivated by the experimental observation of 
Aharonov-Bohm oscillations and persistent currents^ 
In the strictly ID limit, the quantum ring can be stud- 
ied by applying both a discrete model, such as the Hub- 
bard modelj^ and continuum model approaches. In the 
strongly interacting limit, a correspondence is found be- 
tween these two sets of results — the electrons become lo- 
calized and the spin-state of the system becomes that of 
the anti-ferromagnetic Heisenberg Hamiltonian.— £ Elec- 
tron localization causes the system to become a 'rigid- 
rotator'. In this case, only the rotational degree of free- 
dom exists so that a change in the magnetic flux results 
in a change in the angular momentum state of the ring. 

The fractional periodicity seen in the Hubbard 
model and continuum model solutions has been well 
documentedfi^ii^ Physical insight into these proper- 
ties can be obtained from the ground-state energy Bethe 
Ansatz results for the ID Hubbard chain with magnetic 
fluxpi The ground-state energy as a function of flux, 
<&, consists of a sequence of parabolic segments having 
& /N e periodicity. Here, <& = h/e, is the flux quantum 
and N e is the number of electrons. In the strongly inter- 
acting limit, the ground-state solution at a given value of 
the flux can be obtained by creating a single hole in the 
magnon sea. As the flux is increased, the hole is found 
to move from one Fermi point to the next. The U — > oo 
Bethe Ansatz solutions of the Hubbard model show that 
the energy has a (M — N e -^-) 2 flux dependence^ The en- 
ergy minima are found at $/$ = M/N e , where M is the 
angular momentum. The minima of the parabolas and 
fractional periodicity are related to the change in angu- 
lar momentum state as the flux is increased through the 



ring. In the U — » oo limit, the ground-state Bethe Ansatz 
solution becomes that of the anti-ferromagnetic Heisen- 
berg chahiji thereby justifying both the rigid-rotator de- 
scription and anti-ferromagnetism seen in the continuum 
model results^ 

The present work explores the fractional periodicity in 
flux-penetrated t — t' — U ring systems. The t — t' — U lat- 
tice model is a Hubbard model with an additional next- 
nearest-neighbor hopping, t' , termii 2 - For any ID system 
with nearest-neighbor hopping, the ground-state is un- 
magnetized when there is a real and particle-symmetric, 
density-dependent interaction^ Thus, these criteria are 
met for the single-band Hubbard model, but not for the 
t — t' — U system due to the particle-hole symmetry break- 
ing t' termJ^ In the latter case, no definite order is en- 
forced to the particles, and, thus, the possibility of a fer- 
romagnetic state arises. Studies on the ID, t—t' — U chain 
have indeed verified an extensive ferromagnetic phase for 
this system^ Although much work has been done to un- 
derstand the t — t' — U model, little is understood about 
the t — t' — U flux-penetrated ring. Recent investiga- 
tions of the kinetic t — t' (U = 0) system with magnetic- 
flux penetration have indicated an order of magnitude 
increase in the persistent current due to the addition of 
the next-nearest-neighbor hopping term^ To date, no 
investigations on the role of interaction effects in flux- 
penetrated, extended-hopping systems have been made. 

The physics of the flux-penetrated t — t' — U ring will 
be shown to explain the continuum model solutions for 
a circularly symmetric, 2D hard-wall quantum dot. In 
the quantum dot, the strong electron-electron interac- 
tions confine the electrons to a ring-like geometry of fi- 
nite width. The continuum model solution for the flux- 
penetrated system is shown to have ferromagnetic cor- 
relations, with fractional periodicity in the ground-state 
energy as a function of the penetrating flux. The ferro- 
magnetism seen in this system is in contrast to the anti- 
ferromagnetic ground-state result obtained in the limit of 
the purely ID ring. The t — t' — U model, which is con- 
sidered to be the minimal model for ferromagnetism, is 
therefore thought to be an appropriate choice for prob- 
ing the ferromagnetic behaviour seen in the continuum 
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model results. By comparing the lattice model and the 
continuum model solutions, criteria can then be estab- 
lished for obtaining a correspondence between these two 
sets of results. These criteria can then be used to ob- 
tain a microscopic understanding of the ferromagnetism 
in the dot, as well as to explain the fractional periodicty 
in the ground-state energy results. 

The paper is organized as follows. Section II describes 
the lattice model and theoretical method. This is fol- 
lowed by Section III, the results and discussion, which is 
divided into two parts. Part A of Section III details the 
quarter-filled, flux-penetrated t — t' — U ring system. The 
flux-penetrated Hubbard model is first discussed as this 
provides an appropriate context by which to understand 
the flux-penetrated t — t' — U results. Specific criteria 
are then established to explain the onset of fractional 
periodicity seen in the t — t' — U ring as a function of 
t'/t and U/t. Part B of Section III shows an applica- 
tion of the physics of the flux penetrated t — t' — U ring 
by determining the correspondence between it and the 
continuum model results for a few-particle, 2D hard-wall 
quantum dot. The criteria for this correspondence are 
given together with a discussion of the underlying phys- 
ical mechanisms which lead to these results. Section IV 
concludes the main findings of this work. 



II. LATTICE MODEL & THEORETICAL 
METHOD 



The Hamiltonian for the single-band, flux-penetrated 
t — t' — U ring is 

H = - «E( C Hl/»^ + ff. C .) 

icr 

+ ^U l n^n %i . (1) 

i 

Here, c| CT (c^) creates (destroys) an electron with spin 
a = {T, 1} at site i in the system and = c\ a Ci,j is the 
number operator. The parameters t and t' are nearest 
and next-nearest neighbor hopping terms, respectively. 
The L-site ring system is penetrated by magnetic flux $, 
which is encapsulated in the hopping part of the Hamil- 
tonian, such that, 
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L $ ' 



(2) 



The interaction term is the local Coulomb interaction be- 
tween opposite spins, the Hubbard U . All parameters are 
specified in arbitrary units of energy and relative to the 
nearest-neighbor hopping parameter, t. Fig. 1 shows the 
geometry of the t—t' — U ring with flux penetration. Note 
how the t'/t interaction distorts the chain and introduces 
a quasi-2-dimensional aspect to the model. 




FIG. 1: Geometry of the t — t' — U ring with penetrating 
magnetic flux. In accordance with Ref . the combination of 
the nearest and next-nearest neighbor hopping terms is shown 
to distort the ID geometry, resulting in a system consisting 
two coupled ID chains. To illustrate this feature, the two 
rings are displaced along the direction of the magnetic field. 
This is done for visualization purposes only, and does not alter 
the physics of the model. 



The basis functions for the lattice model are defined 
using an occupation number basis, 



\ipa) = \n a i-\...n aN f,n a i l ...n aNl ). 



(3) 



The ground-state eigenstate \^)i a ttic 
combination of these basis functions, 



is then a linear 



HS 

Y&)lattice = 22 a a \^a), 



(4) 



where a a are the coefficients of the basis 
HS = G£ T ) (n, 



states and 

) (jv ) ^ s ^ ne s ^ ze °^ t ne Hilbert space. Here, 
N e -\ and 7V e | refer to the number of spin-up and spin- 
down electrons, respectively. As only small system sizes 
are considered, the eigenproblem remains computation- 
ally tractable and can be solved by exact diagonalization 
using the ARPACK solver The ground-state energy E$ 
and persistent current I are investigated as a function of 
$, where 



W) = 



d£o($) 
9$ ' 



(5) 



To determine the magnetic state of the system, the total 
spin S is calculated. The local moment, (Sf) = 
is used to measure the degree of spin localization. For 
a quarter-filled chain the local moment has a maximum 
value of |. 



III. RESULTS & DISCUSSION 

A. Properties of the Flux-Penetrated t — t' — U 
Chain as a Function of U/t and t'/t 

Before discussing the physics of the flux-penetrated 
t — t' — U ring, it will be necessary to review the flux 
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FIG. 2: (a) Ground-state energy as a function of flux for an 
8-site Hubbard system at quarter-filling (2 f, 2 j), = 20, 
U/t — 200 and 1000. (b) Persistent current as a function of 
flux for the U/t = 1000 result. 



penetrated Hubbard chain as this will provide an impor- 
tant base comparison for the t — t' — U extended model. 
In Fig. 2(a), the ground-state energy as a function of the 
flux is shown for an 8-site Hubbard ring at quarter-filling. 
The results are for three different interaction strengths, 
U/t = 20, 200 and 1000. The periodicity in the ground- 
state energy is found to change from $ to 0> o /N e as 
the interaction strength is increased to large values and 
approaches perfect fractional periodicity as U — ► oo. At 
large interaction strengths, e.g., at U/t = 1000, the angu- 
lar momentum values at the minimum of the parabolas 
follow the expected $/$ — M/N e relation. Thus the 
M = 0, 1, 2, 3 and 4 angular momentum states, for ex- 
ample, correspond to $/$ = 0, 0.25, 0.5, 0.75 and 1.0 
respectively. 

In Fig. 2(6), the persistent current is shown for the 
U/t = 1000 Hubbard ring. The change in positive to neg- 
ative current corresponds to the minima of the parabolas 
in Fig. 2(a). The jump from minimum to maximum cur- 
rent is indicative of the change in angular momentum 
state of this system. The amplitude of the persistent 
current is seen to be almost constant as a function of the 
penetrating flux. This near-perfect fractional periodic- 
ity corresponds to the system being an "effective rigid- 
rotator". At U — oo, perfect fractional periodicity and, 
hence, full rigid rotation would occur. In this case, the 
amplitude of the persistent current would be constant. 

Fig. 3 shows the ground-state energy for the quarter- 
filled t — t' — U system as a function of flux, with U/t = 8 
and t' /t = —0.05. The application of the small t'/t term 
is found to induce marked differences in the ground- 
state energy as a function of the flux strength, thus, 
demonstrating the sensitivity of the solutions to small 
t'/t perturbation. An interesting deviation from the 
Hubbard results occurs — namely, an early onset of the 



FIG. 3: (a) Ground-state energy as a function of flux for 
an 8-site t — t' — U system at quarter-filling (2 j,2 J.), with 
t'/t = —0.05. The values of U/t starting from the bottom 
curve are 0, 1, 2, 4, 6, 8, 10, 20, 40.2 and 1000. U/t = 40.2 is 
specifically shown as this indicates the onset of effective rigid 
rotation, and hence near perfect fractional periodicity, (b) 
Ground-state energy as a function of flux, with focus on the 
U/t = 20, 40.2 and 1000 results. 



fractional periodicity with respect to U/t. This effect is 
clearly evidenced in Fig. 3(6) by the appearance of the 
$/$ = ±0.25 parabolas in the U/t = 20 result. At 
U/t ~ 40.2, the t — t' — U system becomes an effective 
rigid-rotator, which can be seen in the near-perfect frac- 
tional periodicity of the U/t — 40.2 solution. In compar- 
ison, for the Hubbard ring, larger values of the Hubbard 
U (i.e., U/t > 1000) were required in order to achieve this 
effect. The Hubbard results show that the fractional pe- 
riodicity improves as a function of increasing Hubbard U 
strength. In the t — t' — U system, however, the fractional 
periodicity acquires its optimal state at moderate values 
of the Hubbard U. Increasing U/t for the t—t' — U system 
destroys the near-perfect fractional periodicity due to the 
'freezing' of the energetics in the $/3> = ±0.5 result as 
a function of increasing U/t. In this case, the energy cor- 
responding to the ground-state solution at $/$ Q = ±0.5 
remains the same for U/t > 40.2. 

Inspection of the total spin as a function of U/t at 
$/$ = ±0.5 shows that the system has undergone a 
magnetic transition at U/t ~ 40.2, from a S = to 
a 5 = 2 fully-polarized ferromagnetic spin state. The 
ferromagnetic transition can be confirmed by comparing 
the energies corresponding to the lowest energy S = 
and 5 = 2 spin states at $/$ = ±0.5, as a function 
of increasing U/t. At U/t — 20, the ground-state has 
S = with an energy gap of 0.12 to the next state, 
which is 5 = 2. At U/t ~ 40.2, the lowest energy 
5 = and 5 = 2 spin states are near degenerate — the 
5 = 2 ferromagnetic spin state being the ground-state 
solution. Away from this point at U/t = 1000, the fer- 
romagnetic 5 = 2 ground-state persists, with the magni- 
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FIG. 4: (a) The local moment (Sf) as a function of U/t for 
an 8-site t — t' — U system at quarter-filling, (2 f, 2 J,). The 
two solutions are the $/$ = 0.5, t'/t = and t'/t = —0.05 
results, as indicated in the legend. Note the discontinuity 
in the t'/t — —0.05 solution, which occurs at the onset of 
effective rigid rotation (U/t ~ 40.2) and corresponds to the 
5 = 0^2 ferromagnetic transition. 



tude of the energy difference between the lowest energy 
5 = and 5 = 2 spin states, again, being 0.12. The 
freezing of the energetics at = ±0.5 is therefore 

related to the removal of Hubbard U effects in the fully- 
polarized (ferromagnetic) state — in this case, no double 
occupancy and hence Hubbard U contribution can oc- 
cur. The onset of effective rigid rotation and near-perfect 
fractional periodicity corresponds to a system where the 
spins are well localized. Thus, the local moments for the 
5 = 2 ferromagnetic spin-state at $/$ = ±0.5 reach 
the maximum value of (Sf) = 0.375. At U/t = 40.2, 
the local moment results are 0.3748 for $/$ = and 
$/$ = ±0.25. These results are close to the maximum 
value of (Sf) = 0.375. At U/t > 40.2 the system con- 
tinues to evolve as a function of U/t for all states other 
than $/$ = ±0.5. Fig. 4 summarizes the local mo- 
ment results for the t — t' — U system at $/$ = 0.5, 
for the t'/t = (Hubbard) and t'/t = -0.05 cases. The 
ferromagnetic transition for the t'/t — —0.05 system is 
indicated by a change in spin-state to the fully-polarized 
5 = 2 state, together with a corresponding jump in the 
local moment to the maximum value of 0.375. Compared 
to this result, the t'/t = (Hubbard) solutions show a 
smooth evolution of the <&/<& = 0.5 state with respect to 
the local moment solutions and no change in the 5 = 
spin-state of the system. 

Fig. 5 shows how the effective rigid-rotator feature of 
the t — t' — U system evolves as a function of t'/t. In- 
creasing the magnitude of t'/t to 0.1 decreases the ef- 
fective rigid rotation value of U/t to U/t ~ 16.9. The 
onset of effective rigid rotation again corresponds to a 
ferromagnetic transition to the fully-polarized (5 = 2) 
state and freezing of the energetics in the $/$ D = ±0.5 
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FIG. 5: (a) Ground-state energy as a function of flux for an 
8-site t — t' — U system at quarter-filling, (2 |, 2 J.), with 
t'/t — —0.1. The values of U/t starting from the bottom 
curve are 10, 16.9, 20, 30, 40, 50, 100, 200, and 1000. 



TABLE I: Parameter values corresponding to the onset of 
effective rigid rotation, t'/t and U/t, together with AE max 
and the local moment results, (Sf), for the 8-site quarter- 
filled t — t' — U system. The first result is for the Hubbard 
model (t' = 0). For t' 7^ 0, the U/t values also indicate the 
point of transition to the ferromagnetic 5 = 2 state, which 
occurs at &/& = 0.5. 

t'/t ~ U/t AEmax (Sf)$/$ o =0 (5? )$/$„=0.25 {Sj)$/& o =0 .5 
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0.3671 
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0.3750 



results. Consequently, the local moment solutions at 
$/$ D = ±0.5 with U/t — 16.9 have maximum values 
of 0.375, with local moment results for the Q/Q = 
and ±0.25 solutions being 0.3741 in these cases. 

In general, calculations for the t — t' — U system as a 
function of t'/t reveal an inverse relation between t'/t and 
U/t with respect to the onset of near-perfect fractional 
periodicity and effective rigid rotation. Table I shows the 
point of onset of effective rigid rotation as a function of 
t'/t and U/t. Note, this onset also corresponds to the 
ferromagnetic 5 = 0^2 transition in the <&/<!> = 0.5 
result. The inverse relation between the t'/t and U/t 
parameters seen in Table I and corresponding ferromag- 
netic transition are in accordance with the known prop- 
erties of the quarter-filled t — t' — U chain — namely that 
the critcal parameter setting for ferromagnetism occurs 
when U ~ li'l" 1 — see Ref.dl 

The AE max values in Table I are defined as the 
maximum energy difference between the minima of the 
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parabolas at $/$ D = 0, 0.25 and 0.5. At low values of 
t'/t, the corresponding AE max values are small, mean- 
ing that the system has near-perfect fractional period- 
icity. At t'/t — —0.05 and —0.1, for example, the val- 
ues for AE max are 0.001 and 0.003 respectively. The 
relative closeness in energy in these solutions can be ex- 
plained by a minimum Hubbard U contribution to the 
ground-state results occurring at $/$ — 0, 0.25 and 
0.5. Inspection of the highest weighted basis states in 
the ground-state wavefunctions show that at the point of 
effective rigid rotation there are no double occupancies, 
and, hence, no significant Hubbard U contribution to the 
ground-state results. Increasing \t'/t\ adds addtional t'/t 
kinetics thereby causing an increase in AE max . This is 
also reflected in Table I in the corresponding reduction 
in the local moment results. 

In summary, effective rigid rotation in the t — t' — U 
ring can be defined by the following three criteria, 

(1) near-perfect <S> /N e fractional periodicity, with en- 
ergy minima occurring at values of $/$ Q = M/N e . These 
features are characterized by the relative closeness in the 
energy minima of the parabolas as defined by AE max 
and can be accounted for by the freezing out of the U/t 
energetics at $/& = ±0.5, together with a minimum 
Hubbard U contribution at other $/$ = M/N e values. 
(2) An S = — ■> 2 ferromagnetic transition, which occurs 
at $/$ = ±0.5, and (3) maximum or near-maximum 
values in the local moments determined at these points. 

The following trends can also be seen in the persistent 
current results. Comparison of the persistent current for 
the Hubbard (Fig. 6(a)) and t - t' - U (Fig. 6(6)) sys- 
tems show an increase in the number of peaks found in 
the J($) curve. In particular, the U/t = 20, t - t! - U 
solution shows double the number of peaks compared to 
the Hubbard solution. Early onset of effective rigid ro- 
tation is also evident, and this is characterized by equal 
amplitude peaks in the /($) function. Of significance 
are the similarities in the results for the Hubbard system 
at U/t = 1000 (Fig. 6(a)) and the t - t' - U system at 
U/t = 40.2 (Fig. 6(b)) — both showing effective rigid ro- 
tation. Characteristic of the t — t' — U system is that 
it undergoes a ferromagnetic transition leading to the 
freezing of the energetics at $/$ = ±0.5. In addition, 
at <&/$ ^ ±0.5 the energetics continue to evolve beyond 
the point of effective rigid rotation. This leads to a dis- 
tinctive feature in the persistent current curve, namely in 
the formation of signal 'clusters' consisting of the same 
number of peaks as there are particles in the system. 
An example of such a cluster is shown in Fig. 6(b) for 
U/t = 1000 at -0.5 < $/$ < 0.5. 



B. Correspondence Between the Flux Penetrated 
t-t' -U Ring and Few-Particle 2D Hard- Wall 
Quantum Dot 

A circularly symmetric few-particle hard- wall quantum 
dot (QD) is an example of an extended quasi-ring-like 
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FIG. 6: Persistent current as a function of flux, for the 8-site, 
quarter-filled (a) Hubbard chain with U/t = 20 and 1000 and 
(b) t - 1' - U system at t'/t = -0.05 and U/t = 20, 40.2 and 
1000 respectively. Note that the 7 = axis has been shifted 
in order to separate the results for comparison purposes. 

system. Such structures can be formed in semiconduc- 
tor heterostructures by etching techniques or electrical 
gating of the two-dimensional electron gas (2DEG). In a 
typical example of the GaAs/ALjGai-^As heterostruc- 
ture, the trapped electrons can be modelled as two- 
dimensional droplets of charges with renormalized mass, 
to* = 0.067m e , and dielectric constant, e = 12.7eo- 

A conventional choice for a model Hamiltonian for this 
system is 

i— 1 i<j J 

where V(r) is the external confinement potential. In the 
low-density regime and at high magnetic field, the inter- 
actions between the charges are significant. While a real- 
istic interaction potential is much softer due to screening 
effects, a long-range Coulomb interaction whose strength 
is parameterized by the dimensionless scaling factor C 
will be assumed. Additionally, an homogeneous exter- 
nal magnetic field B = V x A, perpendicular to the 
2DEG plane will be included. The confinement poten- 
tial is taken to be of the form 

V(r) = 0, r < R 

= oo, r > R, (7) 

where R = 5a*, is the radius of the dot and a*, is the effec- 
tive Bohr radius. The choice of the hard- wall potential 
is motivated by the aim of having extended ring-shaped 
electronic structures with ferromagnetic ground-state so- 
lutions. This potential has been used in Ref. [l6| to study 
the effect of confinement on the magneto-optical spec- 
trum of a QD, and in Refs. H3 and [la to study the elec- 
tronic structure of a QD in a strong magnetic field, using 
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a mean-field density-functional method. The most com- 
mon choice of potential for quantum dots, the parabolic 
potential, is not a feasible choice as the ground-states 
are not ring-shaped in the regime of interest. In addi- 
tion, the ground-states are not ferromagnetic without a 
strong Zeeman term, despite the effects of Landau level 
mixing 2 ^. This can be seen, e.g., in Ref. [l^, which shows 
exact diagonalization results for electrons in a parabolic 
confinement, treated in the lowest-Landau-level approx- 
imation. The Zeeman coupling is therefore neglected in 
this work to better elucidate the interesting fact that 
the spin-polarized solutions are caused by the electron- 
electron interaction and not by the Zeeman term. 

The ground-state properties of the QD system are de- 
termined numerically by diagonalizing the Hamiltonian 
in Eq. ([6]) in a basis of spin-dependent Fock states^ built 
from a truncated set of lowest-lying single-particle eigen- 
states of the system.— The z-component of the spin, S z , 
and the angular momentum, M, are found to be good 
quantum numbers for this basis. The functional form 
of the energy eigenstates can be obtained analytically, 16 
however, this has explicit dependence on the respective 
energy eigenvalue, which has to be solved numerically. 
The matrix elements of the Coulomb interactions for 
pairs of single-particle states are also evaluated numeri- 
cally. 

Fig. [7] shows the ground-state energy results and the 
radial electron densities obtained from the eigensolutions 
ofEq. ©, with corresponding (M, S) values. InFig.[7Ja) 
the angular momentum of the ground-state energy solu- 
tions can be seen to increase as a function of the increas- 
ing total flux, $ = 7r|B|i? 2 , in the dot. The electron 
density distribution (Fig. [7fb)) assumes a quasi-ring-like 
behavior, with the maximum density occurring near the 
mid-radius. As the flux through the dot and hence the 
angular momentum is increased, the electron distribu- 
tion forms a tighter ring-like structure. For example, at 
M = 10, the radial electron density tends to zero at r = 
and at the edge of the dot at r = 5ao- 

Further evidence of the quasi-ring-like behavior can be 
seen in Fig. which shows the ground-state phase dia- 
gram for the dot as a function of <E> and C. Fractional 
quasi-periodicity in the ground-state parameters occurs 
as a function of the flux for strong values of $ and C. 
Consistent with the flux-penetrated ID Hubbard ring, 
the ground-state total spin varies periodically, with the 
angular momentum increasing in minimal steps of one. 
The spin-polarized ground-states appear in angular mo- 
mentum steps of N e , which is similar to parabolic QD's 
in the regime of small electron numbers. 22 

Interestingly, the total spin oscillation as a func- 
tion of M for the hard-wall QD (Fig. 8) is not the 
same as that reported by Koskinen et al. for the 
quasi-one-dimensional rings^ — the latter system corre- 
sponds to the flux-penetrated anti-ferromagnetic Heisen- 
berg model, whereas the hard-wall QD corresponds to 
the flux-penetrated ferromagnetic Heisenberg model. An 
important question then arises as to what are the micro- 



(a) 




FIG. 7: (a) Ground-state energies and (b) radial electron den- 
sities for the 4-particle hard-wall quantum dot, with M = 
6, . . . , 10. The quantum numbers of the levels are shown in 
(a) as the pair (M,S). The cyclotron energy, hoj c /2, of the 
zero-point motion has been substracted from these graphs. 



scopic mechanisms for these differences between the ID 
ring and the hard-wall 2D systems? 

In order to address this question, the continuum model 
solutions as a function of M and S are compared with 
those from the t—t'—U lattice model, considered to be the 
minimal itinerant model for ferromagnetism in ID sys- 
tems. For convenience, the continuum model results are 
chosen at the flux values shown by the arrows in Fig.[7f a). 
Since the lattice basis is not complete in the continuum 
description, a direct comparison is not possible. Instead, 
the many-body wavefunction ^(xi, . . . , xjv) is projected 
onto a tensor-product basis of delta functions 

\j,a z ) = S(x - x*)^^, (8) 

representative of the lattice basis, with the resulting Fock 
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FIG. 8: Ground-state phase diagram of a 4-particle hard-wall 
quantum dot with R — 5ag as a function of the external mag- 
netic flux, $, and effective interaction strength, C. The shade 
of the filling represents the spin quantum number, from dark 
to light S — 0, 1, 2, and the integers represent the angular mo- 
mentum M of the state. Fractional quasi-periodicity appears 
at strong $ and C. 



expansion normalised to unity. Here, the lattice index j 
runs over the lattice sites, and x* is the real-space co- 
ordinate corresponding to the jth site. The choice of 
x* is not unique, as the only constraint is that the sites 
should form a ring-like lattice. To obtain a unique lattice, 
the radius is chosen to maximize the amplitude of the 
many-body wavefunction with its arguments assuming 
a regular polygonal configuration. Once the wavefunc- 
tion is expressed in the lattice basis, the overlap between 
the continuum and corresponding lattice eigenstatc solu- 
tions, contl^ lattice) , can then be calculated. 

One of the important aspects of this work is to under- 
stand the role of the electronic itineracy with respect to 
the magnetism seen in the continuum model results. To 
determine the relationship between these properties, the 
N e particle continuum model solutions have been pro- 
jected onto a discrete quarter-filled lattice model, having 
(2 x N e ) lattice sites. To justify this choice of lattice, 
a specific example is given of the 4-particle continuum 
model results projected onto an 8-site, quarter-filled lat- 
tice. Here the maximum density droplet (MDD) state is 
considered, which refers to the lowest angular momentum 
state of the QD when the electrons are in the lowest Lan- 
dau levels The angular momentum of the MDD state 
is defined as M M dd = N e (N e - l)/2. For the N e = 4 
particle QD system, M MDD = 6. Inspection of the nor- 
malized ground-state set of eigenvector coefficients deter- 
mined from the continuum model projection onto the dis- 



TABLE II: Overlap results between the 8-site, quarter-filled 
t — t' — U ring and 4-particle continuum model system. The 
lattice model parameters are U/t = 23 and t'/t = —0.17. 



Lcont 


($/$o)cont 


($/$o)iat*ice 


% Overlap 


6 


5.5 


1.5 


98.4 


7 


6.5 


1.75 


98.2 


8 


7.0 


2.0 


97.8 


9 


7.5 


2.25 


97.5 


10 


8.0 


2.5 


97.1 


TABLE III: Overlap results between the 10-site, quarter-filled 


t-t'- 


U ring and 5-particle continuum model system. The 


lattice model parameters are 


U/t = 13 and t'/t 


= -0.19. 


Lcont 


($/$o)«mt 


(<£/$o)la.ttice 


% Overlap 


10 


6.9 


2.0 


98.5 


11 


7.3 


2.2 


98.1 


12 


7.7 


2.4 


97.8 


13 


8.1 


2.6 


97.7 


14 


8.5 


2.8 


97.6 


15 


8.8 


3.0 


97.4 



crete lattice demonstrates that the most highly weighted 
basis states for the -M MDD = 6 ground-state are those 
which are of Heisenberg-type. For the quarter-filled, 8- 
site chain with S z = 0, these highly weighted states 
correspond to those of type (| J, f j 0) and 
(t | I J. 0), with perturbations thereof. There 
are 12 states of this type within the Hilbert space, which 
consists of a total of 784 basis states. For the nor- 
malized wavefunction, the corresponding weightings for 
each of these states is equal to 0.018. The second high- 
est weightings are associated with basis states of type 
(TT | | 00) and (fj. | J. 00), together with per- 
turbations of these, adding up to 96 basis states in total. 
These states have associated weightings of 0.005 each. 
Note specifically that the second highest weighted states 
have nearest-neighbor spins. This justifies the use of the 
extended, quarter-filled lattice and brings to question the 
role of the kinetics with respect to the ferromagnetism 
seen in the continuum model result. 

To determine the correspondence between the t—t' — U 
model and continuum model ground-state wavefunctions, 
the overlaps, con t\^ lattice) , have been calculated. Both 
the 4-particle, S z — 0, and 5-particle, S z = 1/2, con- 
tinuum model cases have been considered. These cases 
were compared to quarter-filled 8- and 10-site t — t' — U 
systems, respectively. The results for the maximum over- 
laps between the continuum and lattice model results are 
shown in Tables II and III. Here, the (3>/$o)cont values 
denote the total amount of magnetic flux that is penetrat- 
ing the dot. For the 4-particle case, these values corre- 
spond to the flux values indicated in Fig. [T^a) . A similar 
choice of flux values was made for the 5-particle system. 
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FIG. 9: Ground-state energy as a function of flux for an 8-site, 
quarter- filled t — t' — U chain with t'/t = —0.17. The values of 
U/t starting from the bottom curve are 3, 5, 7.28, 10, 23, 50, 
100 and 1000 respectively. Note the point of effective rigid 
rotation occurs at ~ U/t = 7.28. 



The maximum overlaps were determined for the first pe- 
riod in M, commencing at the MDD states, Mmdd = 6 
and Mmdd = 10, for the 4- and 5-particle systems, re- 
spectively. The results were obtained by careful scanning 
through the U/t and t'/t parameter ranges. 

Maximum overlaps as a function of M for the 4- 
particle system were found to occur at U/t — 23 and 
t'/t = —0.17 and for the 5-particle system at U/t = 13 
and t'/t = —0.19. The dominant difference in these re- 
sults lies in the strength of the Hubbard U. A larger U/t 
value for the 4-particle system is consistent with what is 
known in the continuum model description of quantum 
rings, namely that these systems become more interac- 
tion dominated when they contain fewer particles. Given 
the similar magnitude of t'/t in both the 4-particle and 
5-particle cases, it seems apparent that the role of this 
parameter is therefore to simulate the correct extended 
geometry of the ring. As M is increased, the contin- 
uum model systems become more ring-like, losing some 
of these extended system qualities. Thus, for larger M, 
the maximum overlaps for both the 4- and 5-particles 
cases are found to slightly decrease. 

In Figs. 9 and 10 the ground-state energy results as a 
function of flux for the 8-site and 10-site t — t' — U ring 
systems are shown with the choice of t'/t values being 
the same as those which were required to achieve maxi- 
mum overlap. The values of U/t required for maximum 
overlap occur after the point of effective rigid rotation in 
both of these cases. From the lattice model perspective, 
this result shows that the system has undergone the fer- 
romagnetic transition required to initiate effective rigid 
rotation. The total spin progression for both the contin- 
uum and lattice model systems as a function of M there- 
fore follows that of the flux-penetrated ferromagnetic 








FIG. 10: Ground-state energy as a function of flux for the 
10-site, quarter-filled t-t' -U chain with t'/t = -0.19. The 
values of U /t starting from the bottom curve are 5, 6.47, 
10, 13, 20, 50, 100 and 1000 respectively. Note the point of 
effective rigid rotation occurs at ~ U/t = 6.47. 



Heisenberg model, namely, (M, 5)4 = (6. 2) — > (7, 1) — > 
(8,0) -> (9,1) -> (10,2) and (M, 5) 5 = (10, f) -» 
(11, |) -> (12, i) -» (13, i) - (14, §) - (15, |), for the 
4-particle and 5-particle cases. The high correspondence 
between the continuum and lattice model systems, being 
between 97 — 98%, demonstrates clearly that, similar to 
the t — t' — U system, the additional degree of kinetic 
freedom and extended ring-like geometry are the impor- 
tant factors which are responsible for the ferromagnetism 
seen in the continuum model results. 

In summary, to encapsulate the correct physics and the 
kinetics in the continuum QD system, the t — t' — U lattice 
needs to be double in size and hence quarter-filled with 
respect to the total particle number. To obtain maximum 
overlap, the lattice model needs to have evolved beyond 
the point of rigid rotation, which is denoted by the on- 
set of a ferromagnetic phase, thus corresponding to the 
ferromagnetism seen in the continuum model results. 



IV. CONCLUSION 

The flux-penetrated t — t' — U, ID system at quarter- 
filling has been investigated as a function of t'/t and U/t 
parameterization. The results indicate an onset of effec- 
tive rigid rotation in this system, which occurs for moder- 
ate U/t values, and coincides with a ferromagnetic tran- 
sition. This model has been used to explain the essential 
physics of the continuum model results for a few parti- 
cle, 2D quantum dot. For 4- and 5-particle cases, the 
maximum overlaps between the continuum and lattice 
model wavefunctions give between 97 — 98% correspon- 
dence. The ferromagnetism, which is seen in the contin- 
uum model solutions, can be explained by the additional 
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kinetics and extended ring-like geometry. The results 
suggest the possibility of anti-ferromagnetic to ferromag- 
netic switching from strictly ID to extended quantum 
ring systems. In addition, a decrease in the order of 
magnitude of the Hilbert space size in the lattice model 
compared to the continuum model indicates that the lat- 
tice model may be a more efficient method of determining 
the ground-state properties and essential physics of ex- 



tended ring QD systems. 
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